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ABSTRACT 

Recent data on the D* masses confirm the QCD prediction based on the very general 
assumption [1] that the spin-dependence of the quark couplings to photons and gluons have 
the same structure and differ only in the value of the coupling constants, This prediction 
with no free parameters and using the H and H* masses as input follows from noting that the 
charmed antiquark c and the strange diquark (ss) have exactly the same couplings to both the 
photon (Q=-2/3) and to the gluon (both are color antitriplets) . New predictions are obtained 
with this approach. Corrections to the naive model are investigated and shown to be small. 







Recent very precise measurements of the isospin splittings in the charmed pseudoscalar 
and vector meson multiplets [2] combined with other experimental masses [3] confirm a nearly 
model-independent QCD prediction [1], 



Sm(D) - 5m(D*) 5m(E) - Sm(E*) , 

10.3 ± 0.6 ± 0.4 x 1(T3 = , - \= \ r^rf = a , J x = 15 ± 4 x 10S 1 

v ; M(D*)-M(D) M(E*)-M(E) K J KJ 

where 5m{h) denotes the mass splitting in a given hadron isospin doublet denoted by h. This 
prediction with no free parameters follows from noting that the charmed antiquark c and the 
strange diquark (ss) have exactly the same couplings to both the photon (Q=-2/3) and to the 
gluon (both are color antitriplets). When a light u or d quark is added to a c to make a D 
or D* or added to a (ss) pair to make a E or S* the strong and electromagnetic contributions 
to the hyperfine splittings have the same ratio in both systems. The hyperfine interaction of 
both of these systems with an additional light quark thus changes by the same ratio when 
the light quark flavor is changed from u to d. This spin dependence of hadron isospin mass 
splittings has been used by Sakharov[4] in the DGG model[5] to obtain an estimate for the 
ratio OLstrongjoL. It has also been used with new data[6] to obtain meson decay constants. 

The above derivation is more general than the previous versionfl] which used additional 
assumptions to treat other hadron masses. Our very general assumption that the spin- 
dependence of the quark couplings to photons and gluons have the same structure and differ 
only in the value of the coupling constants can be applied simply to all color singlet hadrons 
which are bound states of a single u or d quark and some isoscalar combination of quarks or 
antiquarks denoted by X with electric charge Qx and spin Sx 7^ to obtain the relations 



5m(D) - 5m(D*) _ 5m(E) - 5m(E*) _ 6m{E b ) - 5m{E* b ) 
M(D*) - M(D) ~ M{E*)-M{E) ~ M(S*) - M(E b ) ^ ^ 

5m(K) - 5m(K*) 5m{B) - 5m(B*) 5m(E c ) - 5m(E* c ) 



M{K*) - M(K) M(B*) - M(B) M(E* C ) - M(E C 



(26) 



Until recently the errors on the experimental values of the charmed meson masses were 
too great to allow a significant test of the prediction (1). With present data the main source 
of error is seen to arise from E* masses, not the D*'s. Although the errors are still large 
~ 30% the agreement is significant. This is seen by noting the very large difference between 
the ratios (1) and the experimental [3] value —(7 ± 3) x 10^3 for the analogous ratio (2b) 
for kaons. This difference is expected since the electric charges of the strange and charmed 
quarks have opposite signs while the color couplings are the same. Thus the relative signs of 
the electromagnetic and strong contributions to the hyperfine splittings are opposite in the 
two cases. 
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We now examine these relations at a deeper level beginning with the DGG constituent 
quark model [5] in which the hyperfine interaction has the color-spin properties of one-gluon 
exchange, is inversely proportional to the product of quark masses and is treated as a first- 
order perturbation. The most naive version apparently neglects all dependence of hadron 
wave functions on quark masses and assumes that both hyperfine interactions and quark 
magnetic moments scale with the same effective mass parameters. Spectacular success has 
been obtained with this naive approach in using a ratio of roughly 3/2 between effective quark 
masses of strange and nonstrange quarks to explain the observed ratio of hyperfine splittings 
between strange and nonstrange mesons and between strange and nonstrange baryons as well 
as the ratio of A and nucleon magnetic moments. [5,7], 

i k Hp rn s M(A) - M(N) 

1 - 5 = -3^ = ^ = 6 - M(A) + M(iV) +1 = L5 (3a) 
M(A) - M(AQ m. M{p) - M(tt) 

M(E*)-M(S) m u M(K*)-M(K) { ' 

However, the experimental hyperfine splittings in the charmed D mesons do not show this 
factor 3/2. Experiment gives 

Vl± = = 1.01 + 1.5 (3c) 

m„ M(D;) - M(D.) r v ' 



Some insight into this paradox is obtained by noting that increasing a quark mass not only 
decreases the strength of the hyperfine interaction but also increases the value of the wave 
function at the origin and therefore increases the matrix element of the interaction. The two 
effects are in opposite directions and which is dominant is not clear a priori. The result (3c) 
suggests that the two effects may cancel in the charmed case. But the problem remains why 
ignoring wave function effects gives such good results in the lighter quark sector. Furthermore 
this model has as yet no rigorous justification from QCD and the exact meaning of constituent 
quarks and constituent quark masses remain unclear. 

Although detailed models can be constructed to resolve the paradox (3), we consider here 
only the general derivations and implications of the relations (2) which depend on the very 
small u — d mass difference. These relations can be reasonably derived by first-order pertur- 
bation theory and should provide simpler tests of the underlying dynamical assumptions than 
other relations connecting states with larger mass differences. We investigate various wave 
function effects neglected in the simple derivation and note that special problems arise be- 
cause there are two independent perturbations. Both the hyperfine interaction and the mass 
difference are assumed to be small. First-order perturbation theory in a single perturbation 
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automatically takes into account first-order changes in wave functions, and there is danger of 
double counting if one also includes wave function effects in addition to the normal pertur- 
bation results. But with two independent perturbations it is necessary to specify which first 
order wave function effects are automatically included and which are not. 

We first obtain some systematic features of the isospin dependence of hyperfine splittings 
expressible as inequalities confirmed by the systematics of the experimental data. These follow 
from the simple assumption that the effect on the wave function is not strong enough to change 
the sign of the hyperfine splitting. We then show that in a simple constituent quark model 
the relation (1) is independent of the structure of the diquark. Finally we examine the wave 
function effects in a consistent double-perturbation framework, discuss their implications for 
the relations (2), and show that their contribution to (1) is still small, even though they may 
significantly affect relations like (3). 

Using language borrowed from heavy quark symmetry [8], let (Xq)* and (Xq) denote 
the two states with "parallel" and "antiparallel" spins of a light quark u or d denoted by q 
and "isoscalar brown muck" denoted by X. The two states respectively have total angular 
momentum Sx + 1/2 and Sx — 1/2. The mass splitting between such a pair of states is 
assumed in the DGG model [5] to be due to a hyperfine interaction which has contributions 
both from the strong color interaction and from the electromagnetic interaction. We weaken 
the DGG assumptions by not necessarily assuming a constituent quark structure for X nor 
that the strong hyperfine interaction is inversely proportional to quark masses. We only 
assume that the strong hyperfine interaction is greater in the (Xu) configuration than in the 
(Xd) configuration because m u < and similarly for the {Xu)* and (Xd)* configurations. 

It is convenient to define the generalized expression for the ratios appearing in eq. (1) 



The strong hyperfine interaction between X and q in the ground state (Xq) configuration is 
always attractive. It lowers M(Xq) and lowers M(Xu) more than it lowers M(Xd). Thus 
the contribution of the strong hyperfine interaction to the ground state isospin mass splitting 
has the same sign as the mass difference contribution. In the excited (Xq)* configuration 
the strong hyperfine interaction between X and q is always repulsive. It raises M[(Xq)*] 
and raises M[(Xu)*} more than it raises M[(Xd)*]. Thus the contribution of the strong 
hyperfine interaction to the excited state isospin mass splitting has the opposite sign to the 
mass difference contribution. 



F(Qx) 



5m(Xq) - 5m[(Xq)*] 
M[(Xq)*} - M(Xq) 



(4) 



Thus 



S s m(Xq) - 8 s m[(Xq)*} > 



(5a) 



where 5 s m denotes the strong contribution to the isospin splitting. 
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For Qx < and in particular for the case Qx = —2/3 the coulomb and color electric 
interactions in the (Xu) configuration are both attractive; i.e. they both have the same 
sign. Therefore the corresponding magnetic interactions have the same sign. In the (Xd) 
configuration the signs of all electromagnetic interactions are opposite to those in (Xu). Thus 
the electromagnetic hyperfine contributions lower M(Xu) and M[(Xd)*] while raising M(Xd) 
and M[(Xu)*]. Their contributions to the isospin splitting thus also satisfy the relation (5a). 
For Qx > and in particular for the case Qx = +1/3 the electromagnetic contributions 
have the opposite sign from the case Qx = —2/3. Their contributions therefore satisfy an 
inequality opposite to (5a). This can be expressed for the general case as 

-Qx{S e m(Xq) - 5 e m[(Xq)*]} > (5b) 

where 8 e m denotes the electromagnetic hyperfine contribution to the isospin splitting. We 
therefore find that for the case Qx < the total isospin splittings satisfy the inequality 

5m(Xq) - 5m[(Xq)*\ > For Q x < (6a) 

while the ratios F(Qx) defined by eq. (4) satisfy the inequality 

< F(-2/3) > F(+l/3) (66) 

These inequalities are satisfied by experiment. 

We now derive the equalities (2) by introducing the following conventionally accepted 
additional assumptions; 

1. Both the strong and electromagnetic contributions to the hyperfine splittings are given 
by the expectation value of the same operator, not necessarily the value of a wave function 
at the origin. They differ only in the values of the coefficients of this operator. 

2. The strong hyperfine interaction operator factorizes into a factor depending only on 
light quark masses and independent of X and a factor depending upon the properties of X 
and independent of light quark masses. This factorization is found in the DGG model [5]. 
The quark masses used may be constituent quark masses, current quark masses or some kind 
of effective quark masses, all of which satisfy the inequality m u < m<j. 

We consider a "light quark symmetry" for states described as a single light quark moving 
through "heavy brown muck" containing no valence light quarks. We assume that the spin 
and isospin splittings can be treated by perturbation theory and therefore are expressible 
as expectation values of operators describing the perturbation all with the same unknown 
brown muck wave function. We use the Feynman-Hellmann theorem to express the first- 
order dependences of the hyperfine splittings on the mass and the charge of the light quark 



4 



as expectation values of the derivatives of an unknown operator with respect to the mass 
and charge of the light quark. We then find that unknown "brown muck" properties tend to 
cancel in the ratio of the isospin difference between hyperfine splittings to the total hyperfine 
splitting. The result is a universal mass formula for this ratio relating all hadrons in which 
the "brown muck" has the same electric charge. 

Under the above assumptions the hyperfine mass splittings can be written 

A hyp M(Xq) = M[(Xq)*\ - M[(Xq)\ = [S(m q ) - Q x Q q E(m q )} ■ (Xq\ V(X) \Xq) (7) 



where q denotes either u or d, S(m q ) and E(m q ) are constants depending on the value of 
the light quark mass m q and independent of X, Q q denotes the charge of the light quark q 
and the operator V(X) denotes the hyperfine interaction which depends upon the properties 
of X but has its dependence on the flavor of the light quark q factored out and included in 
the coefficients. The constants S(m q ) and E(m q ) and the operator V are defined to make 
S(m q ), E(m q ) and (Xq\ V(X) \Xq) all positive. This leads to the known positive value of 
the hyperfine splitting M[(Xq)*] — M[(Xq)] and is consistent with the observation that the 
color and electromagnetic couplings between X and q have the same sign, attractive, when 
QxQq < 0. Since the electromagnetic interaction is much weaker than the strong interaction; 

i.e. a « CUstrong, 

E(m q ) « S(m q ) (8) 

Then to first order in the small quantities m d — m u and E{m q )/ S(m q ), the Feynman-Hellmann 
theorem gives 

A hyp M(Xd) - A hyp M(Xu) = 



dS(m q ) 
dm„ 



+ £S(m q ) ■ (m d - m u ) + Q x E(m q ) 



{Xq\ V(X) \Xq) 



(9a) 



where the parameter 



log{Xq\V(X)\Xq) 



(96) 



expresses the m q dependence of (Xq\ V(X) \Xq). The need for this correction is discussed 
and justified below. Thus 



A hyp M{Xd) - A hyp M{Xu) _ 1 
A hyp M{Xd) + A hyp M{Xu) 2S(m q ) 



dS(m q ) 



+ ^S(m q ) ■ (m d - m u ) + Q x E(m q ) 



(9c) 

When the wave function correction £ is neglected the right hand side of eq. (9c) is seen to 
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depend upon the constituent X only via the electric charge Qx and therefore has the same 
value for all constituents X having the same electric charge. 

This result (9) can be expressed in a similar form with electromagnetic mass differences 
rather than hyperfine splittings. This gives the function F(Qx) defined by eq(4) as a universal 
function for all isospin differences between hyperfine splittings depending only on Qx- 

F{Qx) = ■ ■ (m d - m u ) + Q x E(m q )) -£-(m d - m u ) (10) 

S(m q ) V om q J 

When the correction term £ is neglected eq. (1) is obtained for the case Qx = —2/3 with 
X = c and X = ss. We can also include X = bs and generalize eq. (1) to (2a). Eq. (2b) is 
obtained for the case Qx = +1/3 with X = s, X = b and X = cs. 

We now examine the validity of the assumption used in obtaining the expression (7) for 
the hyperfine splittings which neglected any effect of the structure of the brown muck wave 
function X. We consider a constituent quark model in which X consists of two constituent 
heavy quarks denoted by x\ and x 2 . Then we can write 

A hyp M(xix 2 q) = M[(xix 2 q)*} - M[( Xl x 2 q)] = A s hyp M(xix 2 q) + A e hyp M( Xl x 2 q) (11) 

where 



^hy P M ( x i x 2q) = s(m q ) ^ (xix 2 q\ g a (q)[g a (l)v( Xl ) + g a {2)v{x 2 )\ \ Xl x 2 q) (12a) 

a=l 



A e hyp M{x 1 x 2 q) = -Q q E(m q ) ■ (xix 2 q\ Q\v(xi) + Q 2 v(x 2 ) \xix 2 q) (126) 

E{m q ) is the same as in eq. (7) and s{m q ) is related to S(m q ) by a color factor evaluated below. 
Both are constants depending on the value of the light quark mass m q and independent of 
x\ and x 2 . The operators v(xi) and v(x 2 ) denote the hyperfine interaction at the constituent 
quark level which again has its dependence on the flavor of the light quark q factored out and 
included in the coefficients. The operators g a {q), fl'a(l) an d <7a(2) denote the eight generators 
of SU(3) co i or acting on the particles q, x\ and x 2 respectively. Eq. (12a) can be simplified by 
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using the SU(3) co i or identity valid for any three-quark color singlet state, 

(xix 2 q\ g a (q) + g a (l) + g a (2) \xiX 2 q) = (13a) 

and noting that 

8 8 8 

CM = J2[g a (q)] 2 = X>(1)] 2 = ][>(2)] 2 (136) 

a=l a—1 a—1 

is the Casimir operator of color SU(3) for a single quark or antiquark state. Then 

8 

A s hyp M(xix 2 q) = -(l/2)s(m q ) ^ ( X 1 X 2V\ 9a{q)[ga{^-) + g a (2)} • [v(x\) +v(x 2 )] \xix 2 q) + A s = 

a=l 

= (1/2) s(m q )Cs(q) (xix 2 q\ v(xi) + v(x 2 ) \xix 2 q) (14a) 

where 

8 

A s = -(l/2)s(m q ) ^2 ( x i x 2q\ g a (q)[ga(l) - g a {2)] ■ [v(xi) - v(x 2 )\ \xix 2 q) = (146) 

A e hyp M(xix 2 q) = -(l/2)Q q E(m q ) ■ [(Qi + Q 2 ) (xix 2 q\ v(xi) + v(x 2 ) \xix 2 q) + A e (15a) 
where 

A e = -(l/2)Q q E(m q )(Qi - Q 2 ) (xix 2 q\ v(xi) - v(x 2 ) \xix 2 q) (156) 

Thus 

A hyp M(xix 2 q) = [S{m q ) - QxQ q E(m q )] ■ (xix 2 q\ V(X) \xix 2 q) + A e (16) 

where 

[S(m q ) = s(m q )C s (q) (17a) 

V(X) = (l/2)[v( Xl )+v(x 2 )] (176) 

The hyperfine mass splittings for the three-constituent-quark system thus has the form of 
the "brown-muck" formulation (7) except for a correction term A e which vanishes identically 
for the case Qi = Q 2 relevant to the expression (1) and can be expected to be small in the 
general case since the operator v(x\) — v (x 2 ) is antisymmetric in the particles 1 and 2 and 
has a vanishing expectation value for a symmetric wave function. 
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We can now also see how the relation (7) can break down for models more complicated 
than simple constituent quark models, where the internal structure of the "brown muck" state 
X effects the strong and electromagnetic couplings differently. The expressions (11-17) can 
be generalized to the case where the "brown muck" state X contains an arbitrary number 
of quarks and antiquarks and the result written in the from (16) with parameters A e and 
A s that vanish in some symmetry limit. However the symmetry limit is no longer relevant. 
There are quarks and antiquarks having different electric charges, and no simple permutation 
symmetry for interchanges of quarks and antiquarks. 

We now show that the effect on the relation (1) of the correction term £ in (10) is small 
in the approximation where (Xq\ V(X) \Xq) depends upon m q only via the reduced mass of 
the (Xq) system. This approximation is exact in a two-body nonrelativistic quark model and 
can be expected to be reasonably good in general. For simplicity we assume a power law 
dependence on the reduced mass. In this case we can write 

f = -nj— ■ log — + = — • 1 — ^ 18a 

Om q \m q mx J m q \ m q + m X J 

Then in the DGG model [5] 

f + «M^)=f-- = — 7— < 186) 

byrriq) \ om q J m q m q m q + mx 

The correction term (18a) is seen to be comparable in magnitude to the dominant term, 
but depending strongly on the ratio — ™ q mx ■ When this ratio is unity; i.e. mx « m q the 

correction term vanishes and all the good results of eq. (3) are obtained. When the ratio 
is small; i.e. mx » m q the correction term tends to cancel the direct term and explain 
the failure of (3) for the charmed sector. We do not attempt here to carry this argument 
beyond hand waving. We consider only the relations (2), where the main correction is seen 
to be independent of mx- We now show that the remaining correction is small. Since the 
expression m q + mx is approximately equal to the mass of the hadron (Xq) , we can express 
the difference in the value of the function F(Qx) (4) for two hadron states denoted by A and 
B having the same value of Qx as 

F(Qx)b - F(Q x )a « n{M M~ M ^ B) ■ - ( 19 ) 

For the case of the D and S hadrons of eq. (1) the correction to the strong contribution is 
~ 8% if we take n=l and m q = 360 MeV and the correction to (1) is ^ 4% if the strong and 
electromagnetic contributions are roughly equal. The effect can be considerably larger in the 
comparison of the K and B systems. It will be interesting to check this experimentally. 
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We now return to the validity of the correction term £ in eqs. (9-10) expressing the m q 
dependence of the matrix element {Xq\ V(X) \Xq). The Feyman-Hellmann theorem tells us 
to ignore wave function effects in first-order perturbation theory because all such effects cancel 
as a result of the variational principle. Therefore the first-order change in (Xq\ V(X) \Xq) 
must be canceled by another contribution. This other contribution arises from the very small 
spin-dependent change in the wave function produced by a change in m q because the wave 
function is an eigenfunction of the total Hamiltonian which includes the hyperfine interac- 
tion. This change in the wave function produces a spin — dependent energy change via the 
spin — independent part of the Hamiltonian; e.g. the kinetic energy. The Feynman-Hellmann 
theorem and the variational principle tell us that effects of this kind must exactly cancel the 
contribution to the energy from the change in (Xq\ V(X) \Xq) to first order. This argument 
suggests that the correction term should not be introduced in eqs. (9-10). 

However, there are really four independent matrix elements and two perturbations: the 
mass difference and the hyperfine interaction. We are calculating a second order cross term 
in a double perturbation series. The Feynman-Hellmann theorem cannot be used twice to 
justify using the same zero-order wave function in the same zero-order Hamiltonian for four 
matrix elements. 

We now clarify this point explicitly and define a consistent perturbation calculation for 
the second order cross term. Consider a general hamiltonian having the form 

H = T(m q )+V + — V hyp (20) 
m,q 

where T{m q ) is the kinetic energy of the light quark but can be generalized to include other 
contributions which depend upon m q , V denotes all contributions to the Hamiltonian which 
do not depend explicitly on m q , e is a small parameter specifying the strength of the hyperfine 
interaction, all the explicit dependence of the hyperfine interaction on m q is in the factor 

and the operator Vhyp contains no explicit dependence on m q . We are interested in the 
dependence of (H) on m q and e to first order in the product e • 5m q , 

Using the Feynman-Hellmann theorem in two ways gives 

= <vw> (21a) 

de m„ 




where the variational principle as expressed by the Feynman-Hellmann theorem allows us 
not to consider the derivatives of the matrix elements. We now obtain the crossed second 
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derivative in two ways by direct differentiation of eqs. (21). Here there is no variational 
principle and the derivatives of the matrix elements must also be considered. 



m 



d 2 (H) _ (V hyp ) d 



q dm q de m q dm q 



+ 7^r ■ ( v hyp) ( 22fl ) 



m M = _M + m i./^) (226) 

q dedm q m q q de \dm q ' 



These results (22a) and (22b) both give the same dominant term — used in eqs. (9) 

when wave function effects are neglected. The result (22a) is just the full eq. (9) and justifies 
the use of the correction term (9b) for the change in the matrix element of the hyperfine 
interaction. Combining eqs. (22) gives 

(Vky P )=m q -.( — ) (23) 



dm q y de \ dm q / 

This condition relating the mass dependence of the hyperfine matrix element and the effect of 
the hyperfine interaction on the kinetic energy seems at first very peculiar. But it must hold 
for any model described by the Hamiltonian (20) since it follows from calculating the same 
second derivative in two ways. It may be useful as a consistency check on detailed model 
calculations. 

We now show how the variational principle relates kinetic and potential energies in the 
manner required by eq. (23) in the case of a nonrelativistic quark model and a short-range 
hyperfine interaction which can be written 

r a r a 

(V hyp ) = J r 2 dr\i,{r)\ 2 v{r) « \^)\ 2 J r 2 drv(r) (24a) 



r r a 
o o 

where ip(0) denotes the wave function at the origin and v(r) is a short range potential which 
vanishes for r > r Q and we neglect the variation in the wave function for r < r . The virial 
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theorem gives 



m /^\=-i-. (U = -JL./ r .^ + ±. r ^\ (25) 
q \ dm q I m q re 2m q \ dr m q dr / 

where T re \ denotes the kinetic energy and /i the reduced mass of the relative motion in the 
(Xq) system. Substituting eqs. (24-25) into the condition (23) gives 

M.3fi.» W -±».(,S ,26) 

1*0(0) \ z om q 2m l q yy 2m q de \ dr / 



This condition must be satisfied in any nonrelativistic quark model with a short range hyper- 
fine interaction and is seen by inspection to be satisfied for the case of a logarithmic potential 

model where r • ^ is a c-number independent of e and ^(O)] 2 varies as // 3 / 2 ). 

We thus conclude that the relation (1) is valid under very general assumptions going be- 
yond the naive DGG model and that it is of interest to improve the precision of the mass 
measurements. The remaining relations (2) of this type are more sensitive to wave func- 
tion effects but should still be approximately valid. The way in which wave function effects 
depending upon reduced mass can explain the paradox (3) is qualitatively indicated in eq. 
(18a). But a more quantitative explanation requires a more detailed model and probably a 
relativistic calculation. 
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